Technical  Research  Report 

Jitter  Analysis  of  CBR  Streams  in  Multimedia  Networks 


by  Jia-Shiang  Jou,  John  S.  Baras 


CSHCN  TR  2003-11 
(ISR  TR  2003-18) 


The  Center  for  Satellite  and  Hybrid  Communication  Networks  is  a  NASA-sponsored  Commercial  Space 
Center  also  supported  by  the  Department  of  Defense  (DOD),  industry,  the  State  of  Maryland,  the  University 
of  Maryland  and  the  Institute  for  Systems  Research.  This  document  is  a  technical  report  in  the  CSHCN 

series  originating  at  the  University  of  Maryland. 


Web  site  http://www.isr.umd.edu/CSHCN/ 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2QQ3  2.  REPORT  TYPE 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Jitter  Analysis  of  CBR  Streams  in  Multimedia  Networks 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Defense  Advanced  Research  Projects  Agency, 3701  North  Fairfax 

Drive, Arlington, VA, 22203-1714 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

see  report 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

unclassified  unclassified  unclassified 

7 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


1 


Jitter  Analysis  of  CBR  Streams  in  Multimedia 

Networks  * 


Jia-Shiang  Jou  and  John  S.  Baras 
Department  of  Electrical  and  Computer  Engineering 
and  the  Institute  for  Systems  Research 
University  of  Maryland,  College  Park,  MD  20742  USA 
{jsjou,  baras}@isr.umd.edu 


Abstract — The  performance  of  real  time  applications 
such  as  video  and  voice  streams  relies  on  packet  delay  jit¬ 
ter.  Large  delay  jitter  causes  buffer  overflow  or  underflow 
at  the  receiver  end  and  the  user  encounters  interrupts. 
The  delay  jitter  is  mainly  due  to  the  perturbation  of 
background  traffic  in  the  bottleneck  router.  Fulton  and  Li 
[12]  provide  an  analytical  approximation  for  the  first-order 
and  second-order  statistics  of  delay  jitter.  However,  their 
analysis  is  based  on  a  Markovian  model  of  the  background 
traffic,  which  is  not  quite  suitable  for  Internet  traffic 
and  requires  lots  of  computational  effort.  We  propose  an 
efficient  method  to  predict  the  jitter  variance  of  a  CBR 
(constant-bit  rate)  connection  based  on  the  wavelet  model 
of  the  background  traffic.  The  wavelet  analysis  extracts  the 
statistical  properties  of  background  traffic  and  the  analysis 
result  can  be  used  to  predict  an  upper  bound  for  the  jitter 
variance  of  the  CBR  connection. 

Index  Terms — Delay  Jitter,  Multifractal,  Wavelet. 

I.  Introduction 

EAL  time  applications  such  as  voice  and  video 
streams  have  become  exceedingly  popular  in  recent 
years.  The  performance  of  real  time  connections  is  very 
sensitive  to  the  connection  quality  such  as  the  packet 
jitter  in  the  network.  The  traffic  behavior  of  a  real  time 
application  is  usually  quite  smooth  and  stable.  After 
entering  the  network,  the  real  time  packet  stream  is 
multiplexed  and  shares  the  link  bandwidth  with  the 
background  traffic.  Hence,  the  packet  jitter  is  mainly 
due  to  the  perturbation  of  the  background  traffic  in  the 
bottleneck  router.  Most  jitter  analysis  methods  [12]  [4] 
[21]  arc  based  on  Markovian  models  of  the  background 
traffic.  However,  recent  studies  [25]  [24]  [20]  [14]  on 
Internet  traffic  have  shown  that  the  aggregate  background 
traffic  driven  by  TCP  is  long  range  dependent  and  self¬ 
similar.  Wavelet  analysis  [29]  [13]  also  demonstrates  that 
background  traffic  is  monofractal  (self-similar)  at  large 

(*)  Research  partially  supported  by  DARPA  contract  No.  N66001- 
00-C-8063. 


time  scales  and  multifractal  at  small  time  scales.  The 
traditional  Markovian  traffic  model  is  unable  to  capture 
this  multifractal  behavior  well  and  is  not  a  proper  model 
for  performance  analysis  [19]  [20].  In  order  to  predict 
the  jitter  more  efficiently  and  accurately,  we  applied 
wavelet  analysis  to  characterize  the  background  traffic 
and  propose  an  upper  bound  for  the  jitter  variance  of  a 
constant-bit  rate  connection. 

The  arrangement  of  this  paper  is  as  follows.  In  the 
next  section,  we  briefly  introduce  the  wavelet  analysis  for 
traffic  and  show  the  multifractal  behavior  of  a  real  traffic 
trace.  In  section  3,  the  background  traffic  is  characterized 
by  the  Logscale  diagram.  An  approximation  of  queue 
length  distribution  is  derived  from  properties  of  wavelets 
and  the  Logscale  diagram.  The  upper  bound  of  CBR 
jitter  variance  is  developed  in  section  4.  Simulation  and 
analysis  results  are  demonstrated  in  section  5.  Section  6 
contains  conclusions  and  suggestions  for  future  work. 


II.  Wavelet  Analysis  of  the  Background 
Traffic 

The  wavelet  technique  is  a  multi-resolution  analysis 
tool  widely  used  in  signal  processing  and  data  analy¬ 
sis  [7]  [5].  It  has  remarkable  advantages  in  analyzing 
stochastic  processes  with  long  range  dependence  [1]  [23] 
[2]  [30]  [16]  [15]  [3].  For  instance,  wavelet  analysis  can 
eliminate  the  effect  of  deterministic  trends  hidden  in  ran¬ 
dom  processes  if  the  wavelet  function  is  chosen  properly. 
Given  the  scaling  function  (f>o  and  the  mother  wavelet  Co, 
the  discrete  wavelet  transform  of  the  continuous  time 
process  X(t)  is  defined  as  follows: 

Definition  (Discrete  Wavelet  Transform)  [23]:  Given 
the  scaling  function  fo  and  the  mother  wavelet  ipo,  the 
approximation  coefficients  aj:k  and  detail  coefficients 
dj,k  of  the  discrete  wavelet  transform  of  the  process  X(t) 
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are  defined  as 


aj,k  J 

roo 

1  X(t)cj)jjk(t)dt 

—  OO 

(1) 

dj,k  '■=  J 

/*oo 

f  X(t)il>jtk(i)dt, 

—  OO 

(2) 

where 

fij,kif)  :  = 

2-j/2fi0(2~jt  -  k) 

(3) 

V :  = 

2~j/2fi0{2  ~jt-k). 

(4) 

It  can  be  shown  that  the  fij  k  and  i/>j  kf°rm  an  orthonor- 

mal  basis.  X(t)  has  the  following  representation 

OO 

X(t)  =  ^2  ao,k(/>o,k(t)  +  ^2^2 

(5) 

k 

For  the  discrete  time 

j= 0  k 

process  Xt,  i  =  0,1,2,... 

,,  the 

discrete  wavelet  transform  can  be  implemented  by  the 
fast  pyramidal  algorithm  [26].  To  understand  the  behav¬ 
ior  of  the  traffic  Xt,  we  arc  more  interested  in  the  detail 
process  of  the  discrete  wavelet  transform  djk-  It  is  well 
known  that  the  Internet  traffic  is  long  range  dependent. 
Many  studies  [8]  [17]  [19]  [6]  have  shown  that  this  long 
range  dependence  property  plays  an  important  role  in 
network  performance. 

Definition  (Long  Range  Dependence)  [18]:  A  sta¬ 
tionary  finite-variance  process  Xi  displays  long  range 
dependence  with  parameter  a  if  its  spectral  density  S(lo) 
satisfies 

S(u)  ~  Cf\uj\~a  as  lo  — >  0,  (6) 

where  0  <  a  <  1  and  Cf  is  a  positive  constant.  It 
also  implies  that  the  autocovariance  function  r[k )  :  = 
E[(Xi  —  EXf)(Xi_ |_fc  —  EXi)\  satisfies 

r(k )  ~  Crka~l  as  k  — >  oo,  (7) 

where  Cr  =  Cf2T(l  —  a)sinfKa/2),  and  T  denotes  the 
Gamma  function. 

The  mother  wavelet  fift)  is  usually  a  bandpass  func¬ 
tion  between  ui\  and  tv 9  in  the  frequency  domain.  Note 
that  the  detail  coefficient  dj  &  is  the  output  process  of  the 
corresponding  bandpass  filter.  The  square  of  the  detail 
process  dj  k  roughly  measures  the  amount  of  energy 
around  the  time  t  =  2j  k  A  and  the  frequency  2_Jteo, 
where  A  is  the  unit  time  interval  and  cco  :=  .l>^2  ■ 

Proposition  [18]:  If  a  stationary  finite-variance  pro¬ 
cess  Xj  has  long  range  dependence  with  parameter  a, 
then  the  corresponding  detail  coefficients  dj  ^  have  the 
following  property: 

log2E[d2. ]  ~  ja  +  log2C(a ,  fi0).  (8) 

Note  that  C(a,  V’o)  is  independent  of  the  variable  j. 


This  property  suggests  that  the  parameter  a  can  be 
estimated  by  the  slope  of  the  log2E[dj .]  v.s.  j  plot.  This 
plot  is  named  the  Log  scale  diagram.  One  advantage  of 
wavelet  analysis  is  that  even  when  the  original  process 
Xi  has  long  range  dependence,  its  wavelet  transform 
djj.  still  has  short  range  dependence  if  the  number  of 
vanishing  moments  N  of  the  mother  wavelet  'iPo(t)  is 
chosen  large  enough  (  N  >  a/2  ). 

Definition  [28]:  The  number  of  vanishing  moments  N 
of  the  mother  wavelet  foil)  is  defined  as: 

ftkMt)  =  0,  k  =  0,1,  2,  ...,N  —  1.  (9) 

Proposition  [9]  [18]:  If  the  number  of  vanishing 
moments  N  >  a/2,  then  dj  k  is  stationary  and  no  longer 
exhibits  long  range  dependence  but  only  short  range 
dependence. 

E[dj)kdf  ifc,]  «  |2 -  2J,A/|a-1-2iV  as  \2>k  -  2^'k'\  -►  00, 

where  j  j'  and  k  k'.  This  implies  the  higher  N,  the 
smaller  the  correlation. 

Figure  1  is  the  Logscale  diagram  of  a  real  traffic  trace. 
Veitch  and  Arby  [27]  developed  an  asymptotically  unbi¬ 
ased  and  efficient  joint  estimator  for  the  parameter  a  and 
C(a,'ipo).  They  also  provided  a  closed-form  expression 
for  the  covariance  matrix  of  the  estimator  and  showed  its 
accuracy.  The  Logscale  diagram  not  only  demonstrates 
the  long  range  dependence  property  of  the  traffic  but  also 
extracts  the  second  order  statistics  at  every  time  scale. 

In  this  paper,  we  characterize  the  traffic  in  terms  of  the 
traffic  mean  rate  and  the  Logscale  diagram.  With  the 
Logscale  diagram,  we  are  able  to  predict  the  overflow 
probability  and  the  jitter  variance  at  the  bottleneck  router. 


slope=0.645  H=0.823 


fine  Octave  j  coarse 

Fig.  1 .  Logscale  diagram  of  a  real  trace 


III.  Traffic  Model  and  Overflow  Probability 

The  wavelet  analysis  has  many  advantages  in 
parameter  estimation  and  traffic  analysis.  The 
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background  traffic  in  this  paper  is  characterized 
by  the  Logseale  diagram  and  mean  rate  m.  An  upper 
bound  of  the  overflow  probability  is  provided  by  using 
the  properties  of  wavelets. 


Lemma  Let  a3j:  and  d3g-  be  the  approximation  coef¬ 
ficients  and  the  detail  coefficients  of  the  Haar  wavelet. 
We  have  the  following  relations: 


ttj,2k  — 
aj,2k+l  = 


®7+l,fc  T  dj+l,k 

"71  ’ 

dj+l,k 


(10) 

(11) 


With  the  uncorrelated  assumption  of  aj  j.  and  djj~  for 
every  j,  we  have 


Var[as\  =  ^K+l]  + 


(12) 


Note  that  the  plot  log 2Var[dj]  v.s  j  is  the  Logseale 
diagram.  On  the  other  hand,  let  Aj  be  the  total  arrival 
bytes  in  the  interval  [0, A).  From  the  definition  of  Haar 
wavelet: 


Aj  =  aj2j/2.  (13) 

Thus,  the  variance  of  workload  Var[Aj\  can  be 
computed  recursively  for  all  j  by  using  the  Logseale 
diagram. 

Lemma  Given  the  Logseale  diagram  log-fA ar{d3\  and 
the  variance  of  Aq,  the  variance  of  A  j,  j  =  1,2, ...  is 

V  ar[A3\  =  2^Var[a3] 

V  ar[a3]  =  2Var[a3-\\  —  Var[dj\. 


Assuming  that  Aj  has  the  Lognormal  distribution  for 
all  j  with  mean  Mj  :=  E[A3\  =  m2J A  and  variance 
Vj  :=  V ar[Aj\ ,  the  probability  density  function  of  the 
Lognormal  distribution  is: 

fAj{x)  ■■=  - ~y=  exp[—  ^ln,T  -jhl.  ];  x  >  o.  (14) 

xojsj  2ir  2cr j 

Since  the  rth  moment  of  the  Lognormal  distribution  has 
a  closed-form: 

r1 2a 2 

EAj  =  exp  (rp3  +  ~^L),  (15) 

the  parameters  pj  and  a3  can  be  easily  calculated  by  the 
following  equations: 

9  +  Vi  x 

=  ln(  M2  )»  <16) 

1  a 

2  ' 


Fig.  2.  The  CCDF  of  workload  Aj  j  =  1,2, ...,  7  and  the  fitted 
Lognormal  distribution 


Fig.  3.  The  CCDF  of  workload  Aj  j  =  8,  2, ...,  14  and  the  fitted 
Lognormal  distribution 


Figure  2  and  Figure  3  are  the  complementary  CDF 
Pr[Aj  >  x]  of  the  real  traffic  trace  and  the  corre¬ 
sponding  Lognormal  distribution  with  the  parameters 
estimated  from  the  Logseale  diagram. 

Consider  a  FIFO  queue  with  an  infinite  buffer  size 
and  the  service  rate  C  (bytes/A).  Assuming  that  the 
distribution  of  Aj  is  known  for  all  j,  Riedi  [22]  proposes 
an  upper  bound  for  the  overflow  probability  P[Q  >  b\. 

Lemma  [22]  Assume  that  Ej  :=  {Aj  <  b  +  C2:> A} 
are  independent  to  each  other  and  the  r.v.  Q  is  the  queue 
length  in  steady  state.  An  upper  bound  of  the  overflow 
probability  of  a  FIFO  queue  is: 


P[Q  >  b\  = 


< 


i  -  P[Q  <  B]  «  l  -  P[nf=1Ej\ 

I< 

l-P[E0\l[P[Ej\Ej_1,...,E0\ 

3= 1 

i-nmi 

3=0 


K 

-nm  <B  +  C2PA],  (18) 

3=0 


Pj  =  In  (Mj) 


(17) 


1 
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where  K  is  the  maximum  octave  and  Ilx  A  is  the 
maximum  time  scale. 

Figure  4  shows  simulation  results  of  the  queue  length 
distribution  and  the  upper  bound.  The  simulation  result 
shows  that  this  upper  bound  provides  a  good  approxi¬ 
mation  of  the  steady  state  queue  length  distribution.  We 
will  apply  this  result  to  predict  the  jitter  variance  of  a 
CBR  traffic  at  the  bottleneck  router. 


0  20  40  60  80  100  120  140  160  180 

Queue  Length  q  (KBytes) 

Fig.  4.  The  queue  length  distribution  with  utilization  p  =  0.4. 


IV.  Upper  Bound  for  Jitter  Variance 

In  this  section,  we  propose  an  efficient  method  to 
predict  the  jitter  variance  of  a  CBR  connection  based 
on  the  Logscale  diagram  of  the  background  traffic  at 
the  bottleneck  router.  Consider  the  following  scenario. 
In  Figure  5,  there  is  one  CBR  process  sharing  the 
bandwidth  and  buffer  with  the  background  traffic  at  a 
FIFO  queue. 


Background  Traffic 


this  paper.  The  packet  size  of  the  CBR  process  is  set  to 
be  small  enough  so  that  its  mean  rate  is  negligible  to  the 
bandwidth  at  the  bottleneck  link.  The  statistical  property 
of  the  background  traffic  is  described  by  its  mean  rate  m 
(bytes/ A)  and  the  Logscale  diagram  Lj  :=  log2  V ar[dj] 
of  the  wavelet  analysis.  The  delay  jitter  of  CBR  is 
defined  as  follows: 

Definition:  Let  the  random  sequence  Ii  be  the  interde¬ 
parture  times  of  the  target  process.  The  jitter  is  defined 
as  the  difference  of  two  consecutive  interdeparture  times: 

Ji  '■=  h+ 1  —  h-  (19) 

We  also  define  An^  as  the  total  arrival  bytes  of  the 
background  traffic  in  the  ith  time  slots.  The  duration  of 
each  time  slot  is  T  :=  2nA  sec. 

Lemma:  Let  the  current  time  be  t  =  TT  and  assume 
that  the  current  queue  length  is  q(t)  >  b  :=  In+1  A(C  — 
m).  The  conditional  variance  of  jitter  is 

Var[J\q(t)  >b\  =  -^Var[(An,)i+1  -  Anj)\.  (20) 

Proof:  Without  loss  of  generality,  let  the  CBR  packets 
arrive  at  times  0,  T,  and  IT,  which  have  queuing  delay 
do,  d\  and  dp,  respectively.  As  shown  in  Figure  6,  the 
total  arrival  bytes  of  the  background  traffic  in  the  ith  time 
slot  (t  E  [IT,  (i  +  1)T))  is  Anj.  Since  the  current  length 
g(0)  >6  is  quite  large,  it  is  reasonable  to  say  that  the 
output  link  is  always  busy  during  the  0th  and  1th  time 
slots.  Moreover,  the  buffer  size  is  infinite  so  that  there 
is  no  packet  loss  event.  The  Lindley  equation: 

q(t)  =  max0<s<t[A(t )  -  A(s)  -  C(t  -  s)],\/t  >  0  (21) 

can  be  simplified  as 

q((i  +  1  )T)  =  q(iT )  +  Anj  -  CT.  (22) 

Hence,  the  packet  delay  d,  =  q(iT)/C  and  the  interde¬ 
parture  time  is 


Fig.  5.  The  target  process  and  the  background  traffic 


Fig.  6.  The  arrival  and  departure  time  of  CBR  traffic 

The  CBR  process  periodically  sends  out  a  small  packet 
every  2"  Asec.  Note  that  the  A  is  the  finest  time  resolu¬ 
tion  in  the  wavelet  analysis  and  we  set  A  =  0.001  sec  in 


Ii  —  T  +  dj+i  —  dj 

q((i  +  l)T)-q(iT) 

c 

A 


The  jitter  variance  under  this  condition  is 


Var[J\q(t)  >  b]  =  Var(Ii+1-Ii) 

Var[(An.i+1-An.i)]  ^ 
C 2  ‘  (  V 

According  to  the  definition  of  wavelet  analysis,  one 
may  easily  obtain  the  value  of  Var[(An^  —  A„x))]  from 
the  Logscale  diagram  Lj 
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for  every  j  and  k. 

Hence,  the  conditional  jitter  variance  of  the  CBR 
process  is 

on+l+Z^+i 

Var[J\q(t)  >b]  = - ^ - •  (26) 

On  the  other  hand,  if  the  current  queue  length  is  small 
( q(t )  <  b),  we  assume  that  there  is  at  least  one  idle  server 
event  happening  in  the  next  two  time  slots.  The  simple 
relations  of  eq.(22)  and  eq.(23)  do  not  hold.  Since  there 
is  at  least  one  idle  event  in  this  period,  the  sequence  of 
interdeparture  times  /,  and  the  sequence  of  packet  delays 
di  can  be  treated  as  uncorrelated  random  sequences 
respectively.  We  have  the  following  approximation: 

V ar[J\q(t)  <  b\  =  Var[Ii+1  -  h\  «  2 Var[Ii\ 

=  2V ar[di+\  —  di]  «  4 Var[dj] 

<  max  V  ar[(Aj  —  C22  A)+]. 

C2  0 <j<n  LV  J 

(27) 


Lemma:  Let  A  be  the  Lognormal  random  variable 
with  parameter  (p,  a)  and  d  >  0  be  any  real  number, 
we  have 


£[(-4  -  o()+] 

£[((.4-rf)+)2] 


a/i+o-2/ 2 


-er/c( 


-  dF(d) 


a2/.l+2cr2 


-erfc( 


In  d  —  p  —  a2  ^ 

~3To  ] 

(28) 

In  d  —  p  —  2a2 


V2> 


a 


-  l£_) 

\/2ct 

+  d2F{d),  (29) 


where  F(d)  :=  Pr[A  >  d\. 

The  probability  of  Pr[q{t)  >  b]  is  based  on  the  pre¬ 
diction  of  the  steady  state  queue  length  distribution.  We 
applied  Riedi’s  approach  [22]  to  calculate  the  overflow 
probability. 

Proposition  Let  the  r.v.  Q  be  the  queue  length  in 
steady  state.  From  eq.  (26)  (27)  and  (18),  there  is  an 
upper  bound  of  the  jitter  variance  of  the  CBR  traffic  at 
the  bottleneck  router: 


Var(J)  <  Var(J\Q  >  b)Pr[Q  >  b] 

+Var(J\Q  <b))Pr[Q  <b].  (30) 


V.  Simulation  results 

The  network  topology  in  our  experiment  is  the  simple 
dumbbell  with  a  single  bottleneck  link.  One  side  of  the 
bottleneck  link  consists  of  800  web  clients,  each  client 
sends  a  web  request  and  has  an  Exponential  think  time 
with  mean  50sec  after  closing  the  current  session.  The 


other  side  has  800  web  servers.  The  server  is  running 
HTTP  1.1  protocol  and  has  a  Pareto  file  size  distribution 
with  parameters  ( K=2 . 3/C bytes,  o-l  .3).  The  propaga¬ 
tion  delay  of  each  server  link  is  uniformly  distributed 
and  the  mean  round-trip  time  is  about  128ms.  Note  that 
the  mean  arrival  rate  of  the  web  (background)  traffic  is 
around  1.2  Mbps. 

A.  Jitter  at  a  FIFO  Queue 

Figure  7  compares  the  predicted  jitter  standard  devia¬ 
tion  with  the  simulation  results.  The  target  CBR  process 
has  fixed  interarrival  times  2nA  and  n  =  3,4,  ..,8. 
The  link  utilization  is  about  0.4,  0.6  and  0.8  with  the 
corresponding  bandwidth  C  =  3.0, 2.0  and  1.5 Mbps 
respectively. 

The  simulation  results  show  that  our  prediction 
method  indeed  provides  a  tight  upper  bound  even  when 
the  CBR  interarrival  time  2 "  A  is  getting  larger. 


Fig.  7.  The  standard  deviation  of  delay  jitter  std(J)  v.s.  n  with  a 
FIFO  queue 


B.  Jitter  at  a  RED  Queue 

We  replace  the  FIFO  queue  by  an  adaptive  RED  queue 
at  the  bottleneck  router.  The  adaptive  RED  queue  [11] 
[10]  will  keep  the  average  queue  length  located  in  a 
desired  region  by  randomly  dropping  the  TCP  packets. 
Since  the  queue  length  is  in  the  desired  region,  the 
link  has  a  100%  utilization  and  no  idle  event  happened. 
Hence,  the  jitter  variance  is  bounded  by  equation  (26). 

Figure  8  shows  that  the  prediction  method  also  pro¬ 
vides  a  tight  bound  for  the  jitter  variance  when  a  different 
queuing  policy  such  as  RED  is  employed. 

VI.  Conclusions 

The  performance  of  video  and  voice  streams  mainly 
depends  on  the  delay  jitter  in  the  network.  Recent  studies 
show  that  Internet  traffic  is  multifractal,  which  can  not  be 
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Fig.  8.  The  standard  deviation  of  delay  jitter  std(J)  v.s.  n  with  a 
RED  queue 

modeled  well  by  traditional  Markovian  models.  Instead 
of  using  a  Markovian  approach,  we  applied  wavelet 
analysis  to  analyze  the  background  traffic  at  the  bot¬ 
tleneck  router.  The  second  order  statistical  property  of 
the  background  traffic  is  characterized  by  the  Logscale 
diagram.  Based  on  properties  of  wavelets  and  some 
reasonable  assumptions,  we  provide  an  efficient  and 
accurate  method  to  predict  the  jitter  variance  of  CBR 
traffic.  Our  simulation  results  also  show  that  this  method 
works  well  with  a  different  queuing  policy  such  as  RED. 
In  future  work,  we  arc  considering  the  effects  of  finite 
buffer  size  in  the  FIFO  queue. 
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